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The dark matter environment of the Abell 901/902 supercluster: a 
weak lensing analysis of the HST STAGES survey 
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ABSTRACT 

We present a high resolution dark matter reconstruction of the z = 0.165 Abell 901/902 
supercluster from a weak lensing analysis of the HST STAGES survey. We detect the four 
main structures of the supercluster at high significance, resolving substructure within and 
between the clusters. We find that the distribution of dark matter is well traced by the cluster 
galaxies, with the brightest cluster galaxies marking out the strongest peaks in the dark matter 
distribution. We also find a significant extension of the dark matter distribution of Abell 901a 
in the direction of an infalling X-ray group Abell 901a. We present mass, mass-to-light and 
mass-to-stellar mass ratio measurements of the structures and substructures that we detect. We 
find no evidence for variation of the mass-to-light and mass-to-stellar mass ratio between the 
different clusters. We compare our space-based lensing analysis with an earlier ground-based 
lensing analysis of the supercluster to demonstrate the importance of space-based imaging for 
future weak lensing dark matter 'observations'. 



1 INTRODUCTION 



Observations and theory both point to the importance of environ- 
ment on the properties of galaxies. Early-type galaxies are typi- 
cally found in more dense regions compared to late-type galax- 
ies ([Dressier 1980), galaxy colo ur and luminosity ar e found to be 
closely related to galaxy density jSlanton et al.ll2005D and the frac- 
tion of star-forming galaxies als o shows a strong sensitivity to the 
densi ty on small < 1 Mpc scales jSalogh et al.l2004l ; lBlanton et al.l 
l2006h . Theoretically there are a number of physical mechanisms 



that could cause these effects in dense environments. These pro- 
cesses can change the star formation history, gas content and/or 
morphology of a galaxy through, for example, ram-pressure strip- 
ping ( Gunn & Gott 1972; Larson et al. 1980; Balogh et al. 2003) 
and/or the tidal effects of nearby galaxies (galaxy harassment, 
iMoore et al.l 19961) and/or the tidal ef fects of the dark matter poten- 
tial ( lBekkilll999l ; lMoore et al.lll998l) . They depend differently on 
cluster gas, galaxy density and the dark matter potential however, 
with ram-pressure stripping dependent on the gas distribution com- 
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pared to tidal effects which are dependent on the overall potential. 
A key difficulty in disentangling these effects observationally, is 
that typically the tidal potential is only constrained in a global sense 
through the measured velocity dispersion of a cluster, or a richness 
or total luminosity estimate. This results in an assumed spherical 
tidal potential model that is smoothed over the small scales that are 
relevant for tidal stripping and harassment studies. 

In this paper we study the complex Abell 901/902 super- 
cluster, hereafter A901/2, in the first of a series of papers from 
the STAGE^I collaboration. From a rich multi-wavelength dataset 
the A901/2 supercluster permits a thorough investigation of the 
relationships between galaxy morphology (fr om Hubble Space 
Telescope (HST ) and ground-based imaging, iGrav et aTl 1200 Sl : 
lLaneetalJ|200% . luminosity, stellar mass and colour (from the 
COMBO- 17 survey wi th 17-band optical imaging, IWolf et al] 
120031: 'Borch et al '}'200&), star formati on rates ( from 24/j m Spitzer 
data. Bell et al. 2007), galaxy density jWolf et al. 2005; G ray et al.l 



2004|). and the hot intra-cluster mediu m (from XMM observations. 



Gilmour et ^1120071 : [Gray et al.ll200Sl . One of the key reasons to 



obtain HST imaging of this supercluster was to construct a high 
resolution, reliable and accurate map of the projected total mass 
density distribution. Using weak gravitational lensing techniques 
we are able to reconstruct the distribution of both dark and lumi- 
nous matter and quantify the significance of the structures that are 
seen, updat i ng the previous ground-based weak lensing analysis of 
iGrav et alj HipO^) . This extra dimension to the multi-wavelength 
view of A901/2 will be a key ingredient in future studies where we 
hope to be able to separate the effects of tidal and gas-dynamical 
influence on galaxy formation and evolution. 

Weak gravitational lensing is now a well established method 
for studying the distribution of dark matter Light from distant 
galaxies is deflected by the gravitational effect of the intervening 
structures, inducing a weakly coherent distortion in the shapes of 
galaxy images. The strength of this lensing effect is directly re- 
lated to the projected mass along the line of sight, and it can there- 
fore b e used to map dark matter in dense regions (see for ex- 
ample iGray et al. I I2OO2I: iGavazzi et alj|2004l: iDietrich et al1l2005l : 
IClowe et al.ll200dlMahdavi et al.ll2007b . 

The first weak lensing analysis of A901/2 by ^ Gray et al] 

1 I2OO2I) used deep ground-based _R-band observations from the 
COMBO- 17 survey jWolf et al .12003 ). This analysis revealed three 
significant peaks in the dark matter distribution at the locations of 
the A901a, A901b and A902 clusters, in addition to a low signif- 
icance south west peak co-incident with a galaxy group, hereafter 
referred to as the SW group. This analysis also showed a filamen- 
tary extension between the A901a and A901b clusters. As this fila- 
ment was loc ated across the CC D chip boundary in the mosaic im- 
age, however, iGrav etaljjiooi) could not rule out the possibility of 
this structure originating from residual uncorrected distortions from 
the point spread f unction (PSF) of the telescope and detector The 
iGrav et alj ( I2OO2I) ground-based analysis also reported a candidate 
giant arc. The STAGES HST imaging can rule out this candidate 
arc as a co-incidental alignment of objects. STAGES does how- 
ever resolve several other cand idate arcs ar ound supercluster galax- 
ies, which will be presented in lAragon-Salamanca et al.; C2008.) and 
lGravetalJ l l2008l) . 

Using the accurate photometric redshift information from the 
A901/2 17-band observations of the COMBO-17 survey, where 



1 'Space Telescope A901/902 Galaxy Evolution Survey' (//Sr GO- 10395, 
PI M. E. Gray), www.nottingham.ac.uk/~ppzmeg/stages 



the photome t ric red shift error ct^ ~ 0.02/(1 + 2) for R < 24, 
iTavlor et al.l ( |2004|) extended the iGrav et al.l ( l2002h analysis, by 
creating a three-dimensional reconstruction of the A901/2 dark 
matter distribution. This analysis revealed a previously unknown 
higher redshift cluster located behind A902 that is at z = 0.46. This 
cluster was named, and hereafter referred to as, CB 1 bv lTavlor et al.l 
(2004). We have updated the redshifts of both CBl and A901/2 in 
this analysis based on an improved photometric redshift catalogue 
and the addition of some spectroscopic redshifts. 

In this analysis we revisit the dark matter distribution in 
A901/2 using deep HST observations. The dominant source of 
noise in the weak lensing analysis of clusters is the Gaussian 
noise introduced from the random intrinsic ellipticities of galax- 
ies. Weak lensing maps of dark matter on small scales therefore 
benefit greatly from the high resolution that HST has to offer. HST 
triples the number density of resolved galaxies from which the lens- 
ing signal can be measured, reducing the intrinsic ellipticity noise 
on small scales. In addition, the high resolution space-based data 
permits higher signal-to-noise shape measurements and a narrower 
PSF, thus implying a more accurate PSF correction. 

This paper is organised as follows. In section|2]we describe the 
weak lensing theory that is related to this analysis, and the maxi- 
mum likelihood method that we use to reconstruct the dark matter 
distribution. We describe the data and weak lensing measurement 
method in section |3] We present our results in section O including 
NFW profile mass measurements in section [43] and the dark mat- 
ter reconstruction and ground-based comparison in section 14.21 A 
first comparison of the dark matter, galaxy light and stellar mass 
distribution is presented in section [43] along with mass, mass-to- 
light and mass-to-stellar mass ratio measurements. A more detailed 
comparison of the mass, gas and galaxies of A901/2 will appear in 
a forthcoming analysis. We investigate the significance of the su- 
percluster substructure that is resolved in our dark matter recon- 
struction in section |4!4l and discuss our findings and conclude in 
section [5] Throughout this paper we assume a ACDM cosmology 
with D.m = 0.3, D.X = 0.7, and Ho = 100hkms"^Mpc"\ All 
magnitudes are given in the Vega system. 



2 METHOD AND THEORY 

Gravitational lensing is sensitive to the projected surface mass den- 
sity along the line of sight E(6'), typically denoted by the conver- 
gence K. In the case of a single lens. 



47rG DiDi, 



(1) 



where Di is the angular diameter distance to the lens. Da is the 
angular diameter distance to the lensed source galaxies, and Dis is 
the angular diameter distance from the lens to the source. 

The coherent distortion, or reduced shear g = gi + ig2, that 
is detected in the images of distant sources allows for the recon- 
struction of the projected intervening matter k. sls g — — it), 
and 

K = +'!/',22) , Jl ^ ^{tp,!! - Tp,22) , 72 = V',12 , (2) 

where 7 is the true shear, 7 = 71 + 172, and ifj^ij is 
the second derivative of the lensing potential (see for example 
iBartelmann & Schneide3l200ll) . 

The strength of all lensing distortions is invariant under the 
transformation k — (1 — A)k + A, where A is a constant 
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jGorenstein et al.ll988h . This is known as the 'mass sheet degener- 
acy' implying that all lensing observations are insensitive to a con- 
stant mass sheet across the field of view (A) in addition to a k de- 
pendent rescaling of the 'original' surface mass density. For wide- 
field images of relatively isolated clusters, where k is weak, one can 
signifi cantly reduce this bias using the C,c statistic of IClowe et al.l 
( Il998r) . The (,c{r) statistic gives a model free estimate of the mass 
enclosed within a radius r and is given by 



Cc(ri) = K(r ^ ri) - K(r-2 ^ r < rmax) : 



(3) 



where r2 is defined to be the radius outside which the cluster den- 
sity is expected to be very low, based on initial mass estimates, 
and Tmax is the field-of-view radius. The second term therefore es- 
sentially measures the constant A. In the case of A901/2 we find 
k(15' «; r < 20') = -0.002 ± 0.007 where r is measured from 
the centre of the STAGES mosaic which is centred on the super- 
cluster. This measure is consistent with what would be expected 
from large-scale structure and the NFW multi-halo model of the 
A901/2 supercluster that we develop in section 143] We therefore 
assume a zero mass sheet degeneracy correction in the analysis that 
follows. 



2.1 Dark Matter reconstruction 

In this paper we use a maximum likelihood method to recon- 
struct the surface mass density k. Starting with a 'best guess' 
lKaiser&Squi^ ( ll993h reconstruction, the lensing potential ^ is 
constructed on a pixelised grid and is allowed to vary to produce 
the minimum difference between the reconstructed and observed 
reduced shear field. The benefit of using this method is that a vary- 
ing noise estimate can be obtained across the whole region enabling 
the significance of each structure in the dark matter map to be ac- 
curately quantified. Furthermore it does not rely on the assumption 
that the observed reduced shear g is approximately equal to the true 
shear 7, which for the A901/2 supercluster would introduce errors 
at the ~ 15% level. We smooth the resulting k maps with a Gaussian 
of smoothing scale 0.75 arcmin, which is equal to ~ 90h~^ kpc at 
the supercluster redshift z — 0.165. This smoothing scale provides 
the best trade-off between high resolution and high signal-to-noise. 

We determine the location of peaks from the local maxima 
and minima in the signal-to-noise weak lensing map. Occasionally 
we find two peaks that are separated by less than half the smooth- 
ing radius. These arise from small noise fluctuations on top of a 
larger fluctuations and in these cases we only count a single peak 
with significance given by the maximal peak within the smooth- 
ing radius. Once peaks are detected in a weak lensing mass map 
their significance has to be compared to what is expected from a 
smoothed random noise map, where a 3a noise peak, for example, 
is much more common t han would naively be expected. As shown 
bv lVan Waerbekd ( l2000h . the statistics of peaks in a smoothed pure 
noise map follow the peak statistics of a t wo-dimensional Gaus- 
sian random field teond & Efstathiodll987h . We use both the peak 
signal-to-noise and the radial peak profile to calculate the global 
probability of a detected dark matter p eak arising from noise using 
Equation (45) of IVan Waerbekd ( I2OO0I) . 



2.2 Model-free Mass measurement 

As our dark matter reconstruction reveals structures that are far 
from the spherically symmetr ic simple isothermal sphere and NFW 
models jNavarro et alj|l997h that are often fit to estimate masses 



from weak lensing measurements (see for examDle [Hoekstral2007h . 
our preferred method to measure mass uses a model-free mass es- 
timate. Following the idea of the (c statistic (Equation[3}, we mea- 
sure the mass of structures within an aperture. For the main struc- 
tures in the supercluster we define apertures by the la and 3a en- 
closed regions in the dark matter signal-to-noise maps. In the cases 
of smaller cluster substructure, where the smoothed structures ap- 
pear to be more spherical, we use circular apertures of radius 0.75 
arcmin to match the smoothing scale used in the dark matter recon- 
struction. The 'aperture' mass is given by 



M= Api^ K{x,y) T.C 

aperture 



(4) 



where Api^ is the projected pixel area at the cluster redshift in 
h~^Mpc'^, {x,y) are pixels enclosed by the chosen aperture and 
Ecrit is the critical surface mass density, given in Equation[T] 



2.3 NFW profile model 

The main drawback of using the model-free mass estimate in Equa- 
tion|4]is the inability to separate mass at different redshifts. This is 
because the dark matter reconstruction k measures the projected 
surface mass density along the line of sight. In the case of the 
A901/2 supercluster there i s a higher redshif t z = 0.46 cluster, 
CBl, that lies behind A902 ^Taylor et al.l2004h such that the model 
free mass estimate for the A902 region gives the combined mass 
of A902 and CBl. To obtain separate mass estimates for the A902 
and CB 1 cluster and to enable a comparison to future analyses of 
numerical simulations, we therefore also present mass estimates for 
the dark matter structures in A901/2 using an NFW halo model. 

The NFW halo model has been shown in numerical simula- 
tions to provide a good fit to the spherically averaged profile of all 
dark matter halos irrespective of their mass ( Navarro et al. 199"^. 
The NFW model for the density profile of a halo at redshift z is 
given by 



p{r) = 



Sc Pc{z) 



(5) 



(r/r,)(l + r-/r,)2 ' 

where Sc is the characteristic density, is the scale radius and 
pc(z) is the critic al density given by 3H(2)V87rG. We follow 
iDolag et al.l ( |2004|) defining the virial radius 7-200 as the radius 
where the mass density of the halo is equal to 200Qrn,{z)pc{z), 
such that the corresponding virial mass M200 is given by 



4i" 3 

M2OO — 200Qm.{z)pc{z) — 7-200 • 

As the mass enclosed within a radius R is given by 



M{r R)=4TTScPc{z)rl 



In 1 



R/rs 



(6) 



(V) 



defining the concentration parameter as c 
teristic halo density 5c is given by 

^ 200 n^(z) 

3 ln(l + c)-c/(l + c) ■ 



1 + R/r. 
r20o/rs, the charaC' 



(8) 



For a given CDM cos mology, the halo mass M200 and con- 
centration c are related (Navarr o et 'Zl ll997l :l Bullock et al. I I2OOII; 
pEke et al. 2001 ; Dolag et al. 200^, where the dependence is calcu- 
lated through fits to numerical simulations. In this paper we use the 
rel ationship be t ween halo mass M200 and concentration c derived 
bv lDolaget"ai] ( l2004l) . 

The expression for the weak lensing shear 7 and convergence 
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K in duced by an NFW dark ma tter halo, given in lSartelmannl il996h 
and lWright & Brainerd Jioool) . depends on the redshift of both the 
lens and source galaxies. In this analysis we have accurate redshifts 
for the majority of the A90 1/902 cluster galaxies but no redshift in- 
formation for ~ 90% of our source galaxies as they are too faint 
to calculate a COM BO-17 photometric red shift. The maximum- 
likelihood method of lSchneider & Rixl ( Il997h was designed to take 
advantage o f such a data set for analysing the galaxy-galaxy lensing 
statistically jKleinheinrich et al.,200 6: Heyma ns et a l. 2006), and it 
is this method that we have adapted for cluster lensing and describe 
below. 

For a model cluster density profile, in the case where all 
galaxy redshifts are known, the weak shear 7 and convergence 
K experienced by each source galaxy can be predicted by sum- 
ming up the shear and convergence contributions from all the 
foreground clusters. In this analysis the redshifts of the source 
galaxies are unknown, and we therefore assign those galaxies a 
magnitude-dependent re dshift probability distr ibution p{z, mag) 
given by Equation 15 of iHevmans et al.l (l2005h updat ed with the 
magnitude-redshift relation of Schrabbac k et al.l where the 

average median redshift Zm is given by 



0.29[mF606w - 22] + 0.31 . 



(9) 



We are then able to calculate the expectation value of the observed 
reduced shear (g) through Monte Carlo integration by drawing a 
source galaxy redshift estimate from the distribution ^(2:, mag), 
V = 1..7Vmc times, where A^'mc = 100 in this analysis. Testing 
larger values for A^mc did not change the result. For each z^ es- 
timate the induced cluster lensing shear g" is calculated with the 
resulting mean reduced shear given by 




Figure 1. The tiling pattern of the STAGES ACS observations. Each group 
of data used to make the seven time-dependent PSF models is shown with 
a different grey-scale. The positions of the four main structures are shown 
with a circle of radius 1 arcmin centred on the brightest cluster galaxy in 
A901a (tile 55), A901b (tile 36), A902 (tile 21) and the SW group (tile 8). 
The position of the infalling X-ray group A901a is also circled (above tile 
55). 



{5> = 



A'"mc 



(10) 



The intrinsic source galaxy ellipticity is then calculated, ~ 
e°^'' — g. The distribution of each component of the observed galaxy 
ellipticity is well described, for the STAGES survey, by a Gaussian 
of width (Je — 0.26. As the induced reduced shear g is relatively 
weak, the probability for observing an intrinsic ellipticity of is 
then given by 

s|2- 

(11) 



1 



■ exp 



2a? 



The best-fit dark matter halo parameters are determined by max- 
imising the likelihood L = 11 [P(e^)i] where the product extends 
over all source galaxies i. 



3 THE STAGES DATA 



The STAGES survey jGrav et alJliooj) spans a quarter square de- 
gree centered on the A901/2 supercluster. Imaged in F606W, using 
the HST Advanced Camera for Surveys (ACS), the 80 orbit mo- 
saic of 80 ACS tiles forms the second largest deep image taken by 
HST. A deta iled accoun t of the STAGES reduction method will be 
presented in iGrav et alj i [200^ . It is ve ry similar to the redu ction 
used for the GEMS s urvey discussed in iHevmans et akl ( |2005|) and 
ICaldwell et aj] ilOOA differing only in the dither and drizzle strat- 
egy. For STAGES, each image consists of four co-added dithered 
images combined with a Gaussian drizzling kernel with a r e sultin g 
0.03 arcsecond pixel scale, as suggested bvl Rhodesetai] j2007h . 
STAGES is complemented by 17-band optical imaging from the 
COMBO-17 survey which provides, for galaxies brighter than R = 



24, accurate photometric redshifts with errors ~ 0.02(1 -I- z), 
spectral energy distribution galaxy classification, and stellar mass 
estimates A/, from low resolution 1 7-band spectra fits to parame- 
terise d star formation history models jBorch et al J 20061 : 1 Wolf et al] 
l2004h . 



3.1 Weak lensing shear measurement 

To measure the reduced weak lensing she ar g, we u se the data 
redu ction steps and method described in iRix et al.] j2004) and 
Hey mans et al. j200S ). The shear measurement aspect is based 
on the iKaiser et al ] 7l995l) method. As we are primarily inter- 
ested in the variation in the dark matter map we have updated 
our shear measurement pipeline to maximise the signal-to-noise 
by including a polynomi al fit to the shear seeing correction P'' 
jLuppino & Kaise 3 ll997h as a function of galaxy size. We also 
include the H oekstra correcti o n to the shear polarisability ten- 
sor detailed in IHevmans et alj 1 I2OO6I) . The accuracy of these up- 
dates has been verified using the publicly av ailable suite of simula- 
tions from the Shea r TEsting Programm^f] jHevmans et akllTood : 
iMassev et al.l |2007|) . The modifications successfully reduced the 
noise on the shear measurement, quantified through the root-mean- 
square variation of the measured ellipticity a^, from = 0.31 to 

a, = 0.26. 

We use the same method as IHevmans et ah H200I to account 
for the time variation of the ACS PSF, namely to divide the data into 
sets imaged in a short period of time and assume that the temporal 



www.physics.ubc.ca/~heymans/step.html 
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variation during that time is minimal. The majority of the A901/2 
field was observed in the space of 20 days, with the remaining 10% 
imaged at a later date over the space of 4 days. Owing to the rel- 
atively low galactic latitude of the A901/2 field and the resulting 
high stellar density of 30 — 40 useful stellar images per ACS im- 
age, we are able to split the data into seven groups to achieve good 
temporal sampling of the PSF distortion. This number was chosen 
to balance between the need to use as many ACS images as possi- 
ble to maximise the signal-to-noise on the average measured stellar 
ellipticity as a function of CCD position, whilst requiring as many 
time bins as possible to minimise the temporal variation of the PSF 
pattern. Figure [T] shows the tiling pattern of the STAGES ACS ob- 
servations denoting each group of data that was used to make the 
seven different PSF models. With this semi-time dependent model 
we find and remove temporal variation during the A901/2 observa- 
tions. Averaged across the ACS field-of-view, this temporal varia- 
tion is at the 1% level on the measured stellar ellipticity. As this 
variation is more than an order of magnitude lower than the weak 
lensing signal from the A901/2 supercluster our semi-time depen- 
dent PSF model is more than sufficient for this analysis. We might 
expect to see low-level systematics for the ACS images whose ob- 
servation date is isolated at the start or end of a data group, affecting 
tiles 21, 33, 36, 43, 44, 46, 47, 57, 69 and 72. Indeed in the analysis 
that follows we find > 3 (t B-modes, an indication of systemat- 
ics dCrittenden et alj|2002h . in tiles 21, 33, 35, 36 and 57. These 
residual systematics will be taken into account by including a con- 
servative systematic error term, base d on the B-mode ampli tude, 
in the analysis that f ollows. Note that lSchrabback et al.l ( |2007|) and 
[Rhodes et af] ( |2007|) present significantly more advanced methods 
to model the temporal variation of the ACS PSF designed for the 
detection of the weaker lensing signal from large-scale structure 
which will be investigated further in a future analysis. 

In the time since the AC S observations of the GEMS survey 
used bv lHevmans et al.l ( l2005h . the charge transfer efficiency (CTE) 
of the ACS has degraded significantly. During the CCD readout, 
as the CTE degrades over time, the amount of charge left behind 
increases. This results in image 'tails' developing along the read- 
out direction, with the most severe effects seen in the furthest ob- 
jects from the readout amplifiers. As the amount of charge left 
behind in each charge transfer is independent of the pixel count, 
CTE impacts on the shapes of fainter objects more significantly 
than brighter objects, and hence this disto rtion is not taken into ac- 
count by the PSF correction. We follow iRhodes et aL |(|2003) by 
using an empirical CTE correction ef for the gi shear compo- 
nent, along the readout direction, where e{™ — A Ay/SN. Ay 
is the distance from the readout amplifier and SN is a signal-to- 
noise estimate that we define as t he ratio of the flux and flux er- 
ror measurements from SExtractor iSertin & AmoutsI 19961) . A is a 
normalisation constant derived to minimise the average measured 
shear (gi), where gi = 2(ei - ef ^^)/Tr(PT), ei is the PSF 
corrected galaxy ellipticity and P~' is the shear polarisabilty tensor 
from Luppino & Kaiser ( 1997). For the faintest galaxies that are 
furthest from the readout amplifier, and hence the most strongly 
affected, ef "^^ = 0.02, but on average ef ~ 0.003 which is 
more than an order of magnitude lower than the weak lensing sig- 
nal from the A901/2 supercluster. We measure the average (gi) be- 
fore applying the CTE correction to be 0.004 and after correction 
{gi) — 0.00001. Any residual CTE distortions that remain after 
the correction are therefore very weak in comparison to the super- 
cluster lensing signal. As the original CTE distortion varies across 
the ACS field of view and hence across the STAGES mosaic, any 



residual CTE distortions would however be included in our B-mode 
analysis and hence the errors in the results that follow. 

3.2 Galaxy selection and redshift estimation 

As we are interested in the dark matter in A901/2 at a redshift of 
z = 0.165, we select galaxies that are likely to be at higher red- 
shifts and thus lensed by the supercluster. As the majority of our 
galaxies are too faint to calculate a COMBO- 17 photometric red- 
shift, the best option is to use a magnitude selection based on the re- 
lationship between median redshift z,n and F606W magnitude de- 
rived in Schrabback et al, ( 2007 ) and given in Equation|9] To ensure 
that the majority of objects have Zs > za90i/2> we select galaxies 
with mpBoew > 23, corresponding to a median redshift Zm > 0.6. 
We also include selection criteria chosen to optimise the accuracy 
and reliability of the weak lensing shear measurement, selecting 
galaxies with S/N > 5, magnitude mpeoew < 27.5, and galaxy 
size rh > 3 pixels. Our resulting weak lensing catalogue includes 
over 60000 objects, or roughly 65 galaxies per square arcmin. The 
average galaxy magnitude of this sample is (mpeoew) ~ 25.7, im- 
plying a median redshift Zm — 1.4. A ssuming a redshift distribu - 
tion given by n{z) (X z^exp(-z^^) teaugh & Efstathioij[T993h . 
we estimate a ~ 3% contamination of our source galaxy catalogue 
from objects that are foreground to the cluster. The dilution of the 
signal by foreground galaxies is therefore well within the statistical 
noise of our analysis. 

To calculate the model-free mass estimate in Equation |4] we 
place the background source galaxy sample at one redshift Zs taken 
to be the median redshift Zm — 1.4 of the sources. We choose 
to use the median, as the high redshift tail of the redshift distribu- 
tion of f aint galaxies is poorly constrained observationally (see for 
example [Benjamin et alj200 7h leading to a potentially biased mea- 
sure of the mean redshift. However note that this choice is fairly 
unimportant as for Zs > 1 and zi = 0.165, the redshift of A901/2, 
the important distance ratio Ds/ Dig in Equation[T]is fairly insen- 
sitive to the value of Zs . For example, a large increase of Zs from 
Zs = 1 to Zs = 1.5 increases Ds/Dis by only ~ 6%. Hence 
for this deep analysis, where the majority of sources have redshifts 
Zs > 1, placing all lensed galaxies at one redshift is a good approx- 
imation. 



4 RESULTS 

In this section we present the results of our weak lensing analysis 
of the A901/2 supercluster including cluster mass estimates, a com- 
parison of three different weak lensing dark matter reconstructions 
and a comparison of the resulting dark matter distribution to the 
distribution of light in the supercluster. 

4.1 Mass estimates for NFW profiles 

We use spherical NFW haloes to model the weak lensing shear 
measured in the A901/2 field. We test two different models using 
the method described in section [23] The 'one halo' model centres 
a single NFW halo on the brightest cluster galaxy (bcgB in each 
cluster at 2 = 0.165. The 'two halo' model places a halo at the 
A901a BCG and the location of the infalling X-ray group A901q!, 

In this analysis we define a BCG to be the brightest cluster galaxy within 
an arcminute of the peak of the cluster's galaxy distribution. 
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Table 1. Mass measurements for the A901/2 supercluster assuming the NFW spherical halo model. The 'one halo' model places a single NFW halo at position 
(RA,Dec) centred on the BCG in each cluster The 'two halo' model places a halo at the A901a BCG and the location of the infalling X-ray group A901a, a 
halo at the A902 BCG and at the background cluster CB 1 BCG. and two halos in the SW group, SWa and SWb. There is no motivation to fit the A901b cluster 
with two haloes and it is therefore only listed in the 'one halo' model upper section of the Table. The NFW 'virial' mass M2oo(h~llOi^M0) corresponds 
to a 'virial' radius r200 (h~ikpc) which has an observed angular scale 6200 (arcmin). For comparison with the 1 arcmin aperture model-free mass estimates 
Map in TableEl M(6» < 1') is the mass of the NFW halo enclosed by a 1 arcmin aperture, centred on (RA,Dec). The reduced of the fit is given in the final 
column. 
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Figure 2. The tangential reduced shear distortion as a function of the dis- 
tance from the BCG in each cluster The dashed line in each panel shows the 
best-fitting model profile assuming single NFW dark matter haloes centred 
on each BCG. The solid line shows the best-fitting model profile assuming 
multiple NFW daik matter haloes. The upper panels show the profile ex- 
pected from three NFW haloes centred on the BCG in A901a, the BCG in 
A901b and the X-ray infalling group A901o. A901a and A901o increase 
the large scale single in A901b and vice versa. The lower left panel shows 
the profile expected from two NFW haloes centred on the BCG in A902 
and the BCG of the background cluster CB 1 . The lower right panel shows 
the profile expected from two NFW haloes SWa and SWb in the SW group. 
The model halo parameters are given in Table [T] 



a halo at the A902 BCG and at the background cluster CB 1 BCG 
at z = 0.46, and two halos in the SW group, SWa and SWb. The 
positions of the two SW halos are motivated by the dark matter 
reconstruction presented in section l42l Table [T] lists the resulting 
constraints on the NFW 'virial' mass M200 of each halo, the NFW 
'virial' radius r2oo and the corresponding angular projection of this 
radius on the sky 0200- We find results that are fully consistent with 
the single and mu ltiple simple isothermal sphere halo analysis of 
lTavloretal.l ( l2004.) . For comparison with the model-free mass esti- 
mates in section l431 we also calculate the NFW halo mass enclosed 
by an aperture of 1 arcmin M {6 < 1') using Equation|7] 

Figure[2]compares the measured tangential reduced shear dis- 
tortion gt around each cluster with the prediction from our 'one 
halo' (dashed) and 'two halo' (solid) model, assuming all lensed 
source galaxies are at a single redshift Zs — 1.4. The rotated shear 
gr is found to be consistent with zero on all scales as expected. 
N ote that the model param eters are constrained using the method 
of ISchneider & Rixl jl997h as described in section [23] not as a fit 
to the azimuthally averaged shear gt presented in this Figure. We 
find that both models fit this data equally well but, in the cases of 
A901a and the SW group, NFW haloes are in general a poor fit 
to the data (the reduced xt of ths fit is given in the final column 
of Table [TJ. From this we conclude that for these unvirialised sys- 
te ms, the spherica l NFW model provides a poor fit. Note that using 
the lBullock et al.l ( l200lh relationship between virial mass and con- 
centration results in even poorer fits to the data. Figure |2] is also 
instructive to see the impact of the individual halos on each other. 
The effect of A901a on the profile of A901b (and vice versa) can be 
seen from the increasing signal on large scales in the upper panels 
of Figure|2] The decrease in signal on small scales in the SW group 
data favours the 'two halo' model over the 'one halo' model. 

Comparing the A902 'one halo' and 'two halo' models in the 
lower left panel of Figure [2] shows that the addition of the back- 
ground cluster CB 1 has only a weak effect on the profile of A902. 
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We find that the virial mass of A902 decreases by 11% when the 
CBl halo is included in the analysis (see Table [T](. As CBl and 
A902 are separated by 1.4 arcmin, the mass enclosed by a 1 ar- 
cmin aperture M{9 < 1') is even less effected and decreases by 
6% when the CBl halo is included in the analysis. We therefore 
conclude from this NFW analysis that in the model-free mass es- 
timates that follow, the contribution from CBl to the A902 mass 
cannot be more than a 10% effect which is within our B-mode 
systmatic errors in the analysis that follows. 

4.2 Dark Matter maps 

Figure [5] shows the STAGES dark matter reconstruction of the 
A901/2 supercluster. The upper left panel shows the maximum like- 
lihood reconstruction that clearly reveals the four main superclus- 
ter structures; A901a, A901b, A902, and the SW group. The con- 
tours on this signal-to-noise dark matter map correspond to -4a,- 
2a (dashed), 2a, 4a, and 6a detection regions (solid). Assuming 
constant noise across the image, the scale bar shows the measured 
convergence k. This is a very good assumption except for the edges 
of the map where the noise increases rapidly. 

The dark matter map can be compared to the B-mode or 'sys- 
tematics map' in the lower left panel of Figure [3] Thi s is created 
by rotating the galaxies by 45° jCrittenden et alj 20021) and recon- 



structing the map. As weak lensing produces curl-free or E-mode 
distortions, a significant detection of a curl or B-mode signal in- 
dicates that ellipticity correlations exist from residual systematics. 
Comparing the B-mode map with the contours from the dark mat- 
ter map therefore allows one to assess the reliability of each de- 
tected structure. For the maps shown in Figure [S] a 3a detection 
has K ~ 0.07, although the true significance of any peak in the 
distribution has to be deter mined by comparison to the statistics 
of a random Gaussian field dVan Waerbek3l2000h . For a field this 
size, with the same number of galaxies, ellipticity distribution and 
smoothing scale, smoothed Gaussian noise would produce 2 ± 3 
random > 3a E and B-mode peaks, and Oil random > 3.5a 
E and B-mode peaks, which we discuss further in section l4!4l All 
but one of the three most significant > 3.5cr B-mode peaks can be 
linked to regions where the simple semi-time-dependent PSF mod- 
eling used in this analysis fails, as discussed in section [bTI 

For compariso n with previous analyses , the upper right panel 
of Figure[3]shows a lKaiser & Sauiresl ( ll993 l KS93) reconstruction. 
The main difference seen between our preferred maximum likeli- 
hood reconstruction (upper left panel) and the KS93 reconstruction 
(upper right panel) is the strength and significance of the peaks. In 
the cluster cores, where k > 0.12, the reduced shear g is more than 
15% larger than the true shear 7. Hence the KS93 k reconstruction 
where the reduced shear g is assumed to be equal to the true shear 
7, results in an overestimation of k in these regions. When the maps 
are smoothed, these strong overestimated k peaks are smeared out, 
and in the case of the neighbouring A901a and A901b and a large 
smoothing scale, this smearing could lead to false filamentary fea- 
tures that we start to see between A901a and A901b in the KS93 
reconstruction. It is this effect, in addition to the possible presence 
of residual PSF systematics, that we conclude are responsible for 
the filamentary extension between A901a and A901b seen in the 
COMBO- 17 dark matter reconstruction o f iGravetalJ 12002*). 

The lower right panel of Figure[3]compares our STAGES HST 
dark matter map (contours), with the COMBO- 17 E SQ 2.2m Wide- 
Field Imager dark matter map of iGrav et alj ( I2OO2I) . This compar- 
ison shows excellent agreement in the locations of the dark mat- 
ter peaks. The ground-based map is shown on the same k scale as 



the other maps, but smoothed with a 1 arcmin Gaussian, compared 
to the 0.75 arcmin Gaussian used in the STAGES. Increasing the 
resolution of the ground-based map by narrowing the smoothing 
scale increases the noise in the map, thereby further lowering the 
sigrrificance of th e detected peaks. In the ground-based analysis of 
iGrav et all f20o3) 21 galaxies per square arcmin were used for the 
dark matter reconstruction with a root-mean-square variation of the 
measured ellipticity of a^ = 0.44. This can be compared to the 65 
galaxies per square arcmin used in this analysis with a^ = 0.26. 
The increased number density of objects in this analysis results 
from the high HST image resolution. The reduction in a^ results 
from both the higher average signal-to-noise imaging of the source 
galaxy sample and the higher average galaxy-to-PSF size ratio that 
can be achieved with space-based imaging. For a 1' x 1' area of 
sky, the random intrinsic ellipticity noise on the measured shear is 
cr^(l') = 0.1 for the ground-based analysis and a-,{l') — 0.03 for 
this space-based analysis. As the weak lensing signal that is typi- 
cally detected around z ~ 0.2 clusters is 7 ~ 0.1 and of the order 
of the ground-based noise, this comparison clearly demonstrates 
the need for space-based observations for better than arcminute res- 
olution 'imaging' of the dark matter. 



4.3 A comparison of mass and light 

The distribution of dark matter in A901/2 is found to be very well 
traced by the distribution of galaxies associated with the superclus- 
ter, as shown by Figure|4] This Figure shows the total r-band lumi- 
nosity of cluster galaxies, smoothed on the same scale as the dark 
matter map (shown with contours). Cluster galaxies are identified 
from the ground-bas ed multi-colour CO MBO-17 data using the se- 
lection criteria from I Wolf et al.l ( |2005|) : their photometric redshift 
must lie in the range 0.155 < Zp < 0.185 and their absolute V- 
band magnitude Mv < —17, which corresponds to an apparent 
_R-band magnitude R < 21.5. These criteria were chosen to keep 
field contamination low and cluster completeness high, with asam- 
ple that is 68% complete at this luminosity limit. I Wolf et aljf2005l) 
estimate the level of field galaxy contamination to be 3% for the 
red-sequence ga laxies, and 15% for the blue-cloud galaxies, (see 
IWolf etalj2 005. for more details). 

Immediately in Figure |4] we can see that the most massive re- 
gions are also the most luminous. We also start to see cluster sub- 
structures repeated in both the dark matter and light maps. A com- 
parison of mass and stellar mass in the cluster is nearly identical to 
Figure |4] implying that the stellar mass is also a good tracer of the 
underlying dark matter distribution. Figure|5]shows an 8 x 8 arcmin 
close-up of the four main structures of the A901/2 supercluster. 
This Figure compares the distribution of dark matter (shown con- 
toured) to the lumi nosity weig hted distribution of old red sequence 
galaxies defined in I Wolf eFal.! ( |2005|) . The locations of the BCGs 
are shown with filled diamonds. For A901a and A901b, the max- 
imal peak in the dark matter distribution is practically co-incident 
with the location of the BCG, (within 0.25 arcmin). For A902 we 
find two peaks in the dark matter distribution matching the two 
BCGs. The dark matter peaks are slightly offset from the BCGs 
(0.5 and 1 arcmin) due to the presence of CBl, the background 
cluster at a redshift of 2; = 0.46 whose BCG location is shown 
in the A902 lower left panel of Figure |5] with a star. The cluster 
CBl fills ~ 1 arcmin aperture around the BCG. The NFW 'two 
halo' A902 and CBl model detailed in section 143] predicts a shift 
in the observed A902 dark matter peak by ~ 0.3 arcmin which is 
consistent with what we find in the dark matter map. 

For the SW group, there is again good agreement with the po- 
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Figure 3. The dark matter reconsttiuction of the A901/2 supercluster. The maximum likelihood dark matter map is shown in the upper left panel clearly 
revealing the four main superclusters structures; A901a, A901b, A902 and the SW group. This reconstruction can be c ompared to a|Ka iser & SquiresI j 19931 
KS93) reconstruction (upper right panel), the associated B-mode or systematics reconstruction (lower left panel) and the lGray et al.l )2002l) ground-based dark 
matter reconstruction (lower right panel). The contours enclose the -4(t,-2(t (dashed), 2a, 4cr and 6(7 detection regions (solid) of the STAGES maximum 
likelihood reconstruction (upper right and repeated on the lower left panel), the Bmode reconstruction (lower right panel) and the KS93 reconstruction (upper 
left panel). The scale bar shows the amplitude of the mapped convergence k which has been smoothed on 0.75 arcminute scales corresponding to ~ 90h~^ 
kpc at the supercluster redshift z = 0.165. The edges of the maximum likelihood dark matter map and the corresponding B-mode (left hand panels) go to 
zero as a result of the data mask that is shown in the tiling pattern of Figure[T] 



sition of the peak in the mass distribution and the BCG, although 
for this group there are two local maxima in the dark matter dis- 
tribution. Interestingly there are two distinct groups in the galaxy 
population of the SW group. There is an old red galaxy population 
that surrounds the BCG, as shown in the lower right panel of Fig- 
ur e[5] In a ddition there is a dusty red galaxy population (described 
by I Wolf e t al. 2005, but not shown in the Figure) that exists to the 
east of the BCG and co-incides with the most massive eastern dark 
matter peak (denoted SWb in Table [T]!. A more detailed analysis 
of the interesting relationship between the dark matter environment 
and the different galaxy populations will be presented in a future 



paper. The lower right panel of Figure |5] also shows one case of a 
significant density of old red galaxies without a peak in the dark 
matter distribution. Towards the edge of the STAGES imaging, the 
noise in our dark matter map grows rapidly, and at the location of 
this galaxy group the noise is twice the noise level at the SW group. 
As this group is likely to be less massive than the SW group, which 
is detected at 5(7, we are not surprised that this group is undetected 
in our dark matter map. 

The A901a upper left panel of Figure |5] shows a significant 
extension of the dark matter distributi on, in the directio n of the 
in-falling X-ray group A9Qia found bv lGrav et alj ( 1200^ (shown 
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Figure 4. A comparison of mass and light in A901/2. The mass distribution 
from Figure[3](shown contoured) is compared to the smoothed light distri- 
bution of the cluster galaxies. The scale bar shows the r-band luminosity 
per square arcmin in units of h~^ lO^^Lr©- 





Figure 6. A comparison of mass and the mass-to-stellar mass ratio A//A/» 
in the main structures of the A901/2 supercluster; A901a (upper left), 
A901b (upper right), A902 (lower left) and the SW group (lower right). 
Ic to 7a contours of the signal-to-noise dark matter map shown in Figure|3] 
are drawn over a smoothed mass-to-stellar mass ratio map. The locations 
of the brightest cluster galaxies are shown (filled diamonds), in addition to 
the location of the in-falling X-ray group A901o (filled triangle), and the 
location of the higher redshift z = 0.46 cluster CBl (star). The scale bar 
shows log(Af/Af,). 






Figure 5. A comparison of mass and light in the main structures of the 
A901/2 supercluster; A901a (upper left), A901b (upper right), A902 (lower 
left) and the SW group (lower right). Icr to la contours of the signal-to- 
noise dark matter map shown in Figure [5] are drawn over a smoothed lu- 
minosity map of the old red sequence supercluster galaxies. The locations 
of the brightest cluster galaxies are shown (filled diamonds), in addition to 
the location of the in-falling X-ray group A901a (filled triangle), and the 
location of the higher redshift z = 0.46 cluster CBl (star). Local maxima 
in the dark matter map are shown with a cross. 



with a filled triangle). The peak along the extension, with Kpoak > 
4(7 (shown with a cross), is co-incident with the brightest galaxy in 
the the X-ray group. 

Comparing the local maxima in the A901b and A902 distri- 
bution with the light maps we find that the substructure in the dark 
matter maps are often associated with substructures in the galaxy 
distribution. The only striking discrepancy is a luminous peak to the 
north west of A902, seen in Figure |4] This luminous peak results 
from a single, very luminous dusty red galaxy that is brighter th an 
the BCG and is likely to be infalling on A902 jWolf 



s bngnter th a 
et al.il2005l) . 



In Table [2] we list mass and mass-to-light ratios for the main 
structures shown in Figure |5] As discussed in section 14.11 these 
structures are far from the spherically symmetric NFW models that 
are often used to constrain models. We therefore use a model-free 
mass estimate given by Equation |4l defining the enclosed region 
using the Icr and 3a contours sh own in Figure [5] F or comparison 
with the ground-based analysis of I Gray et alj j2002l) and the NFW 
analysis of section|4jT]we also list the mass enclosed by a 1 arcmin 
circular aperture (denoted 'ap') centered on each clusters BCG. To 
estimate the contribution of systematic error to our mass estimate 
we follow the conservative prescription that is often used in the 
analysis of weak len sing by large-scale structure (see for example 
iBeniamin et al.ll2007b . calculating errors by adding the random er- 
ror (listed as the first mass error in Table [2l( in quadrature with the 
B-mode signal, shown in the lower left panel of Figure[3]and listed 
as the second mass error in Table|2] The systematic error dominates 
the random error in this analysis. 

We find A901a and A901b to be the most massive systems 
in the supercluster with masses ~ 6.5h"HO^^M0 and mass-to- 
light ratios ~ 150 for the full extended Icr region. A902 and the 
SW group have similar masses, roughly half the mass of the A901 
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Table 2. Model-free mass measurements for the main structures in the 
A901/2 supercluster. The enclosed masses M, mass-to-light ratios M/L 
and mass-to-stellar mass ratios M/M^ are given for the regions defined by 
the la and 3cr contours in the signal-to-noise maps. The area of these re- 
gions is given in the second c olumn (square a rcmin). For comparison with 
the ground-based analysis of iGrav et al.l b002D we also list the mass and 
mass-to-light enclosed by a 1 arcmin circular aperture (denoted 'ap') cen- 
tered on the cluster BCG. The quoted mass errors are listed with the ran- 
dom noise error followed by a conservative estimate of the systematic error 
The mass-to-light ratio error and mass-to-stellar mass ratio error include 
only the eiTor on the mass with the random and systematic parts added in 
quadrature. 



but the reader should note the significance of the mass detected in 
these regions (shown contoured) which is less than 2<j in all cases. 

Comparing our convergence «; mass reconstruction results to 
our NFW shear analysis, we find very good agreement between 
the mass measured within 1 arcmin of each clusters BCG which 
provides an important verification of our dark matter reconstruc- 
tion method. Indeed we find this good agreement between the two 
methods continues out to a radius of 4 arcmin, after which contri- 
bution from to the k map from neighbouring clusters complicates 
the comparison. The difference that is seen between the NFW virial 
masses quoted in Table[T]and the masses quoted in Table[2]is only a 
result of the different physical scales probed in both Tables, which 
can be seen by comparing the observed NFW virial scale 6200 with 
the region area quoted in the first column of Table[2l 

In comparing our results to the previous ground-based lensing 
analysis of the A9 0I/2 supe rcluster we must first consider the as- 
sumption made bv lOrav et al. (2002) that k <C 1 and hence 7 = <?, 
where g is the measured reduced shear given above Equation [2] 
For the main structures in A90I/2 this would result in an overes- 
timate of cluster mass by ~ 15%. Taking this overestimate into 
account, our mass estimates are consistent with lGrav et al.l ^20021) 
as can be seen from the space/ground mass map comparison in the 
lower right panel of Figure[3] The mass-to-light ratio measurements 
for A901b and A902 disagree however at the 3(t level. The dif- 
ference arises from improvements in the selection of the cluster 
galaxies from the COMBO- 17 data, in comparison to the previ- 
ous two-band optical cluster selection of Grav et al. (2002). This 
improvement removes the striking differe nce between the c luster 
mass-to-light ratio measurements found bv lGrav et al.l j2003) . Our 
results show a mass-to-light ratio within an aperture of 1 arcmin of 
M{< 1') ~ lOOh^^MoL© to be a good description for all the 
main structures in the supercluster. We find a similar result for the 
mass-to-stellar mass ratio where M/M*(< 1') ~ 25 for all the 
main structures in the supercluster. 



pair at ~ 3.5h"H0^^MQ . We find A901b to be the most extended 
structure in the system, and the SW group is the most compact. 

The mass-to-stellar mass ratios Af/A/, of each structure are 
given in the final column of Tabled These mass rati os were calcu- 
lated w ith a Hubble parameter h = 0.7, assumin g a lKroupa et al.l 
( I1993I) initial mass function i lBorchetal.ll200^ . The results are 
equivalent to w ithi n 10% of the same result derived using a 
IChabriej ( l2003h or a lKroupj ( 1200 ih initial mass function. We find 
mass-to-stellar mass ratios Af/Af, that are simil ar to the ratios 
found for massive ellipticri ggJaxies at this redshift ( Hoekstra et al.l 
I2OO5I : iMandelbaum et alj|2006l : iHevmans et ai]|200^ . although a 
direct comparison is hard to draw as the results from the massive 
elliptical galaxies measure NFW virial mass to stellar mass ratios 
instead of the model free mass ratio estimates that we present here. 

Figure[6]shows the variation of the mass-to-stellar mass ratio 
across each of the main structures in A90I/2 on a log scale, com- 
pared to the mass distribution (shown contoured). Note that nega- 
tive regions in the mass reconstruction have been set to zero in this 
Figure. Moving out from the central BCG (shown with diamonds), 
we find that the mass-to-stellar mass ratio initially increases, as the 
stellar mass decreases more rapidly than the halo mass. Continuing 
out further, the mass-to-stellar mass ratio then rapidly decreases as 
the dark matter mass tends to zero. This Figure shows some regions 
of very high mass-to-stellar mass ratio regions (log A//Af, > 2), 



4.4 Supercluster substructure 

In this section we investigate the lower significance peaks in the 
dark matter distribution that are not associated with the cores of the 
supercluster structures discussed above. Table[3]lists the number of 
local maxima and minima in the dark matter reconstruction for dif- 
ferent significance levels and compares them to what we find in our 
B-mode reconstruction and what we would expect from a smoothed 
random Gaussian field using Equation (41) from Van Waerbek3 
( I2OOOI) . The high significance peaks ACpcak > 4(t are all associ- 
ated with the cores of the four main supercluster structures. How- 
ever we can see that we have a significant number of Kpcak > 2(t 
peaks that cannot be explained by random noise alone. There are 
a comparable number of |Kpoak| > 2cr peaks in the B-mode map, 
but comparing the location of E and B-mode peaks allows one to 
assess the reliability of the lower significance E-mode detections. 

In order to distinguish between noise peaks and true peaks in 
the mass distribution, it is useful to add morphological information 
about the profile of the peaks. The mean p rofile and dispe r sion o f 
a noise peak is given by Equation (47) of IVan WaerbeS ( l2000h . 
Comparing the measured profile around each detected peak with 
the mean noise profile allows for the calculation of the probability 
that a peak with a given signi ficance and shape is a noise fluctu- 
ation, (using Equation (45) of I Van Waerbekd j2000l) ). In Tabled 
we list the number of peaks that have a less than 33% probability 
of being a random noise fluctuation. The result is consistent with 
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'=pcak< '^pcak> «pcak N^..^ (M«0.75')> N^f^^ {M / L) 

-5a -4a -2a 2a Aa 5a la (h-ilOl^Mg) (hMe/Lj.©) 

Noise 0±0 19 ±5 19 ±5 0±0 2<t-3(j 14 0.35 ± 0.04 10 293 ± 24 

3(T-4(T 2 0.57 ±0.06 1 317 ±114 

Signal 1 36 29 7 3 1 

Bmode 16 21 



Signal (p) 1 19 25 7 3 1 

Bmode (p) 10 14 

Signal (b) 1 19 23 7 3 1 



Table 4. Mass measurements for the halo substructure sample, assuming all 
peaks in the matter distribution are associated with the cluster. Over half of 
the peaks have associated cluster galaxies (^''fj^^]^). for which we measure 
a mass-to-light ratio ({M/ L)). 

with cluster light and are unlikely to be caused by noise or system- 
atics, we find two peaks that are also co-incident with small galaxy 
groups of 3-4 galaxies at the same redshi ft z ~ 0.45 . As th is is 
same redshift of the CB 1 cluster found in Taylor et al. we 
are potentially seeing extended large-scale structure at higher red- 
shift which is supported by the fin dings of an optical cluster search 
of COMBO- 17 data in this field jpalter et al.ll2008h . This will be 
investigated further in a forthcoming three-dimensional analysis. 



Table 3. Peak statistics; comparing the number of peaks as a function of sig- 
nificance in a smoothed random Gaussian field (Noise), the reconstructed 
dark matter map (Signal) and the reconstructed B-mode map. We use the 
peak profile to further distinguish between a noise peak and a true peak. 
Signal (p) and B-mode (p) list the number of peaks in our substructure sam- 
ple, where the peak has less than 33% chance of being a noise peak. The 
Signal and B-mode peaks are very rarely co-incident, but as a conserva- 
tive measure we define a high confidence sample, (Signal (b)), where the 
B-mode at the peak must be less than half the amplitude of the signal. 



the difference between the total number of detected peaks and the 
expected number of random noise peaks. 

To define a low-significance 2a" < Kpeak < 4(t substructure 
sample we use high confidence selection criteria where the peak 
must have less than 33% probability of being a random noise fluc- 
tuation, and the B-mode at the location of the peak must be less than 
half the amplitude of the E-mode. The last row of Table[3]lists the 
number of peaks that meet these criteria for different significance 
levels, leaving 16 'substructure' peaks with 2a < Kpeak < 4a. 
Note that the 7 peaks with Kpcak > 4a are all associated with the 
central regions of the four main structures in the supercluster, dis- 
cussed in section|4l as shown by the marked crosses in Figure|5]that 
are enclosed by the 4a contour. 

Assuming all 16 substructure peaks in the dark matter map are 
associated with the supercluster, we can calculate a mass for these 
halos using Equation|4] We find an average mass of M{< 0.75') = 
0.35 ± 0.04h"^10"Mo for the 2a < Kpcak < 3a group, and 
M(< 0.75') = 0.57±0.06h-H0^^Mo for the 3a < Atp^ak < 4a 
group. Table |4] lists the number of peaks that are associated with 
cluster galaxies where L(< 0.75') > lO^'h^^Lr©. We find that 
over half of the peaks are associated with galaxies in the cluster, 
and provide average mass-to-light ratios for these peaks in Table|4] 
The mass-to-light ratio of our 'luminous' substructures are of the 
same order of magnitude as the mass-to-light ratios of the main 
supercluster structures listed in Table |2] 

It is likely that many of the peaks that are not associated with 
cluster light are actually at a different redshift, as the dark mat- 
ter map shows the projected surface mass density along the line 
of sight. We have found one particularly interesting 3.5a peak in 
the dark matter distribution to the south west of A901a, that is not 
co-incident with any cluster light. This peak has a 0.1% chance of 
being a noise fluctuation and is not co-incident with any signifi- 
cant B-modes. Our hypothesis is that this peak is due to a mass 
concentration at a higher redshift than the cluster, supported by the 
presence of a small group of 5 galaxies found within a 0.8 arcmin 
aperture, centred on the 3.5a dark matter peak, which have pho- 
tometric redshifts z = 0.44 ± 0.04. Intriguingly, out of the four 
less significant 2a — 3a dark matter peaks that are not associated 



5 DISCUSSION AND CONCLUSION 

From a weak lensing analysis of deep Hubble Space Telescope 
data, we have reconstructed a high-resolution map of the dark mat- 
ter distribution in the Abell 901/902 supercluster. We find that 
the maximal peaks in the dark matter distribution are very well 
matched to the locations of the brightest cluster galaxies in the most 
massive structures in the supercluster. These structures are A901a, 
A901b, A902 and the South West group, all of which are detected 
in our dark matter map at high significance. 

Owing to the high number density of resolved objects in the 
HST data, we have been able to produce a map with sub-arcminute 
resolution. This has allowed us to resolve the morphology of the 
dark matter structures, finding profiles that are far from the spher- 
ically symmetric NFW models that are typically used to model 
such systems. We find local maxima in the dark matter distribution 
around the main structures, that are also seen in the distribution 
of galaxies. Furthermore we see a significant extension in the dark 
matter distribution around A901a, in the dir ection of an in-falling 
X-ray group called A901a l lOrav et al.ll2008l) . 

We have presented mass, mass-to-light and mass-to-stellar 
mass ratio estimates for each of the main structures, finding A901a 
and A901b to be the most massive clusters in the system with M{< 
1') ~ 2 X 10^^h~^Mo. Contrary to the analysis of lOravet all 
( I2OO2I) we find no evidence for the variation of the mass-to-light 
ratio or the mass-to-stellar mass ratio between the different clusters 
measured in a 1 arcmin aperture 120h~^ kpc) centred on each 
cluster. We have shown the variation of the mass-to-stellar mass 
ratio across the clusters, finding an initial rise in M/Af, as a func- 
tion of distance from the clusters central BCG, followed by a steep 
decrease. 

We have investigated the less significant substructures in the 
dark matter map that are detected at < 4a. Comparing the profile 
of these peaks with what is expected from a random noise peak 
we have selected a sample of substructures where the likelihood of 
those peaks being noise or a result of an imperfect PSF correction 
is low. We find that over half of these peaks are associated with 
galaxies in the cluster, yielding mass-to-light ratios that are com- 
parable to the mass-to-light ratios found in the main structures in 
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the supercluster. The remaining peaks in the distributi on are hkely 
to be associated with galaxy groups at higher redshift jpalter et"al] 
l2008 h. supported by the discovery of several co-incident groups of 
galaxies at 2 '-^ 0.45. 

One interesting r esult of lGrav et alj ( I2OO2I) was a tentative de- 
tection of a filamentary extension between A901a and A901b. We 
do not recover this signal in this analysis however and conclude that 
this feature was a result of residual PSF systematics and the KS93 
mass reconstruction method used in the Gray et al. ( 2002) ground- 
ba sed analysis. A simil ar non-detection and conclusion was drawn 
bv lOavazzi et al.l j2004h on a re-analysis of the tentative lensing de- 
tection of filamenta ry structure in the MS0302-I-17 supercluster by 
lKaiseretal.Ul998l) . These two null results do not mean, however, 
that filamentary extensions of dark matter do no t exist between 
clusters. Instead, as shown by iDolag et alj ^20061) . we are finding 
that intra-cluster filaments are very difficult to detect through weak 
lensing. From numerical simulations, Dolag et al. (2006) determine 
an expected filamentary shear signal from a supercluster filament 
of (? ~ 0.01, which is a factor of three smaller than the noise on 1 
arcmin scales in this HST analysis. To detect a signal of this magni- 
tude would require significantly deeper space-based observations. 
An alternative, that we are currently investigating, is the detection 
of weak gravitational flexion, a third order weak lensing effect that 
will be very effective at probing the sub-structures t hat were re- 
solve d in this weak shear analysis (see for example iBacon et al.l 
l2006l) . and is also a potential way to recover more information 
about intra-cluster filaments. 

The dark matter map presented in this paper will form the ba- 
sis of future studies of galaxy morphology and galaxy type in an 
over-dense dark matter environment. Comparing the results of this 
analysis with the previous ground-based analysis clearly demon- 
strates the importance of space-based observations for future high 
resolution weak lensing dark matter 'observations' of dense envi- 
ronments. 
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